Suppression of decoherence via strong intra-environmental coupling 
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We examine the effects of intra-environmental coupling on decoherence by constructing a low 
temperature spin-spin-bath model of an atomic impurity in a Debye crystal. The impurity interacts 
with phonons of the crystal through Jahn- Teller vibronic coupling. Anharmonic intra-environmental 
vibrational coupling is incorporated through anti-ferromagnetic spin-spin interactions. The reduced 
density matrix of the central spin representing the impurity is calculated by dynamically integrat- 
ing the full Schrodinger equation for the spin-spin-bath model for different thermally weighted 
eigenstates of the spin-bath. Exact numerical results show that increasing the intra-environmental 
coupling results in suppression of decoherence. This effect could play an important role in the 
construction of solid state quantum devices such as quantum computers. 
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Intra-environmental coupling has customarily been ne- 
glected in theoretical models of subsystem-environment 
interaction. The popular spin-boson model |l[ |§1 , for ex- 
ample, assumes that the environment consists of a set of 
non-interacting harmonic oscillators linearly coupled to a 
central spin. The neglect of intra-environmental coupling 
is motivated more by mathematical convenience than by 
physical insight. 

One failing of such models is that intra-environmental 
energy transfer can only proceed by using the subsystem 
as an intermediate. In addition, it is well known that 
the statistical properties of the energy eigenfunctions and 
eigenspectra of strongly coupled (irregular) systems is 
qualitatively different from that of uncoupled (regular) 
systems. The Wigner functions of eigenvectors for irreg- 
ular systems are almost uniform over the energetically 
allowed classical phase space 1^. Those of regular sys- 
tems are more lumpy with energy localized in a subset of 
the available modes |Q. Similarly, irregular eigenspectra 
show level repulsion while regular spectra show level clus- 
tering These spectral signatures have important dy- 
namical consequences For these reasons coupled en- 
vironments may have decoherence properties which differ 
substantially from those predicted by uncoupled oscilla- 
tor models. 

Proposed new technologies such as quantum comput- 
ing , laser control of chemical reactions ||^ , and molec- 
ular electronics all require a qualitative understanding 
of the effects of decoherence and dissipation for experi- 
mental implementation. Predictive theoretical studies for 
such systems would also greatly benefit from dynamical 
methods which accurately include the effects of intra- 
environmental coupling and environmental memory ef- 
fects. 

Unfortunately, exact theories such as the Feynman- 
Vernon influence functional method jlOj and the 
Nakajima-Zwanzig master equation pd] | cannot be eas- 
ily applied. Approximate theories such as the Red- 
field p2[, completely-positive-dynamical-semigroup flS] 



and SRA |l4| master equations need testing against ex- 
act results before they can be applied with confidence. 
Thus, exact numerical calculations for subsystems inter- 
acting with environments of a small number of degrees of 
freedom provide the only practical method for exploring 
the effects of intra-environmental coupling. Such stud- 
ies would also allow us to test the accuracy of existing 
master equations. 

In this manuscript we report exact numerical results 
for the decoherence of a central spin interacting with a 
spin-bath with intra-environmental coupling. The model 
was constructed to represent an impurity in a thermal 
crystalline solid; the impurity being the small quantum 
system and the solid the environment. The first two 
electronic states of the impurity are vibronically cou- 
pled to a number rig of coupled phonon modes of the 
crystal. At low temperature the phonon modes can be 
represented as spin-1/2 modes [|l^ (representing the first 
two states of an oscillator) with frequencies sampled from 
a Debye spectrum with a frequency cutoff (set here at 
lod ~ 1). With anharmonic phonon-phonon coupling 
effects included, our model Hamiltonian then takes the 
form 
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where we arbitrarily chose loq = .8288 as the frequency of 
the impurity, /3 = .01 is the coefficient of a small anhar- 
monic correction, and Aq = 1 and A are the subsystem- 
environment and intra-environmental coupling constants. 
Terms one, two and four of ([l]) represent the uncoupled 
modes of the subsystem (labeled with superscript 0) and 
environment (labeled with superscripts 1 through ris). 
The third term in (|l|) couples the subsystem and envi- 
ronment, while the last term couples the environment 
with itself. Here the sigmas represent Pauli matrices. In 



2 



our units h = 1. [Note that the same Hamiltonian could 
also represent interacting nuclear or electronic spins in a 
solid.] 

We calculated the reduced density matrix p{t) of the 
impurity via the formula 
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trn and |m)e are the energies and eigenvectors of the iso- 
lated environment, and kT is the temperature in units of 
energy. The notation Tre{A} indicates a trace of the op- 
erator A over the environmental degrees of freedom. The 
states |'0jn(i)) are evolved via the Schrodinger equation 
from initial states 
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under Hamiltonian The basis of eigenstates of the az 
operators was used to represent all states. The states |0) 
and |1) represent down and up z-components of the spin, 
respectively. Thus, the subsystem state |1) in Eq. (||) 
means that the impurity is initially in its excited state. 

Equations ^ and (||) represent an impurity in a ther- 
mal solid which is excited by a fast laser pulse just prior 
to time t = Q (i.e., a sudden approximation) which then 
evolves while interacting with phonons in the solid. 

The calculations reported here are for = 12 bath 
spins. The ARPACK linear algebra software |]l6| was 
used to calculate the lowest n^ig = 20 energies and 
eigenvectors of the isolated environment. A tempera- 
ture of kT — .02 was chosen such that no states with 
energy quantum number m higher than Ueig would be 
populated at equilibrium. The numerical solution of 
the Schrodinger ordinary differential equation for |V'm(0) 
was calculated using an eighth order Runge-Kutta rou- 
tine Operations of the Hamiltonian (|l|) on the 
wavevector were calculated via repeated application of 
Pauli matrix multiplication routines. For example 
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for all sets of j; = 0, 1, I 

ji = and ji = if = 1. Thus, an operation of ax 
simply rearranges the components of j?/;). States of the 
basis can be represented by integers j — ji ^- ~\~ ■ ■ ■ ^ 



. + j„^2"= and since integers are represented 



in binary form on a computer, the mapping j — > j' — 
ji +j22 + ...+ ii2^-i + _ + j„,2"=-i under ai*^ can 
be calculated very simply using Fortran binary-operation 



intrinsic functions. Operations for ay 
straightforward. 
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Figure 1; Subsystem entropy S versus time t plotted for dif- 
ferent values of A 
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In Fig. |l| we show the calculated subsystem entropy 
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where det[p(t)] = pii(t)poo(i) - Pw{t)poiit), for various 
values of the intra-environmental coupling A. For A = 
(solid curve) the entropy rapidly approaches its maxi- 
mum value of ln(2) ~ 0.693147. . . which corresponds to 
the diagonalised reduced density matrix 
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As A increases through 1 (long dashed curve), 2 (short 
dashed curve), 4 (dotted curve), and 8 (dot-dashed 
curve), the entropy approaches a smaller asymptotic 
value. This strongly suggests that increased intra- 
environmental coupling suppresses decoherence. 

To confirm this we compare the dynamics of the three 
subsystem spin components 
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calculated for Hamiltonian (|^) with those of the impurity 
evolving in isolation (i.e., Aq = 0). In Fig. ||we show the 
time evolution of the central spin in the absence of inter- 
action with the bath, i.e., for Aq = (solid line) and the 
evolution of the same spin for Ao = 1 and A = (dashed 
line), i.e., when the central spin is coupled to a bath of 
non-interacting spins. As can be seen from Fig. |2(c)| , the 
system undergoes rapid decoherence when the interac- 
tion with the bath is turned on. In Fig. ^ we show again 
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Figure 2: a, y, and z component of the spin versus f for A = Figure 3: x, y, and z component of the spin versus f for A = 8 
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the behaviour of X{t), Y{t), and Z{t) for Aq = (sohd 
hne, isolated spin) and for Aq = 1 and A = 8 (dashed 
line, central spin coupled to a bath of strongly interact- 
ing modes). Here the agreement between the coupled 
and isolated dynamics is much better than in the previ- 
ous case; in fact, the only significant discrepancy appears 
in the behaviour of Z{t), whereas the plots of X{t) and 
Y{t) for Ao = and Aq = 1 are almost indistinguishable. 

Thus, at least in this simple model, increasing anhar- 
monic intra-environmental coupling incrementally sup- 
presses decoherence. There is a relatively simple explana- 
tion for the observed behaviour. Define an environmental 
super-spin with components 

Us rts 
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in terms of which Hamiltonian (0) can be rewritten as 

H = ^ai°)-KM°)+AoafS. 

where Q = X^j^i^j/^-s average environment fre- 

quency and Vj = ujj — f2. When A ^ £7, the last term 
in (^) dominates the bath Hamiltonian. This entails that 
the bath eigenstates are nearly eigenvectors of Y,x ; specif- 
ically, the bath eigenstates with lowest energy correspond 
to the lowest eigenvalues of Tix- Since the initial condi- 
tions (H) are low-energy bath eigenstates, the evolving 
dynamical states \ij}rn{t)) remain close to eigenstates of 
Tlx with small eigenvalues; as a consequence the central 
spin becomes nearly decoupled from the environment and 
its evolution is determined by the effective Hamiltonian 

^ (c.o/2)af + + Aos,)4°) 

where Sx is a (typically small) eigenvalue of T,x ■ 

Our study strongly suggests that intra-environmental 
coupling has as important an effect on decoherence as 
temperature or subsystem-environmental coupling. En- 
vironments of N (coupled oscillator) phonons should 
have similar decoherence properties. Since the Wigner 
functions of the energy eigenstates of strongly coupled 
systems are nearly uniform over the classical energy sur- 
face [|, energy is distributed equally among all modes. 
Thus displacements from equilibrium of any phonon 
mode must be small at low temperature, especially in the 
thermodynamic limit. Since the coupling of an impurity 
to a phonon is through its displacement coordinate, this 
coupling will also be small. By contrast, energy distribu- 
tion in uncoupled oscillator systems is non-uniform and 
may be localized in a small number of modes [Q . Hence 
vibronic coupling to these modes will be strong as will 
decoherence. 



We have obtained similar results for both larger and 
smaller numbers of bath spins and also in the case of 
ferromagnetic intra-environmental coupling. We believe 
that our model is a reasonably accurate representation 
of an impurity in a low-temperature crystal. If such 
strong intra-environmental coupling exists in nature, it 
could be exploited as a platform for quantum comput- 
ing. We are currently developing methods for the study of 
higher-temperature systems where environmental modes 
are modeled by coupled harmonic oscillators. 
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